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Abstract 

Resonant relaxation (RR) of orbital angular momenta occurs near massive black holes (MBHs) where the 
potential is spherical and stellar orbits are nearly Keplerian and so do not precess significantly. The resulting 
coherent torques efficiently change the magnitude of the angular momenta and rotate the orbital inclination 
in all directions. As a result, many of the tightly bound stars very near the MBH are rapidly destroyed by 
falling into the MBH on low-angular momentum orbits, while the orbits of the remaining stars are efficiently 
randomized. We solve numerically the Fokker-Planck equation in energy for the steady state distribution of a 
single mass population with a RR sink term. We find that the steady state current of stars, which sustains the 
accelerated drainage close to the MBH, can be < 10 larger than that due to non-coherent 2-body relaxation 
alone. RR mostly affects tightly bound stars, and so it increases only moderately the total tidal disruption rate, 
which is dominated by stars originating from less bound orbits farther away. We show that the event rate of 
gravitational wave (GW) emission from inspiraling stars, originating much closer to the MBH, is dominated 
by RR dynamics. The GW event rate depends on the uncertain efficiency of RR. The efficiency indicated by 
the few available simulations implies rates < 10 times higher than those predicted by 2-body relaxation, which 
would improve the prospects of detecting such events by future GW detectors, such as LISA. However, a higher, 
but still plausible RR efficiency can lead to the drainage of all tightly bound stars and strong suppression of GW 
events from inspiraling stars. We apply our results to the Galactic MBH, and show that the observed dynamical 
properties of stars there are consistent with RR. 

Subject headings: black hole physics — Galaxy: center — stellar dynamics — gravitational waves 



1. INTRODUCTION 

Galactic nuclei with massive black holes (MBHs) are stel- 
lar systems with relaxation times that are often shorter than 
the age of the Universe. In that case the distribution function 
(DF) of the system may approach a steady state. This steady 
state is determined by the boundary conditions (an inner sink, 
such as the tidal radius of the MBH, and an outer source at 
the interface with the host galaxy) and by the mutual interac- 
tions between the stars themselves. The nature of the "micro- 
scopic" interactions between the stars determines the rate at 
which the system relaxes to its steady state. 

With the notable exception o f A^-body simu lations (e.g. 
Baumgardt, Makino & Ebisuzaki 2004al| 2004bl Preto, Mer- 
ritt & Spurzem r2004t Merritt & Szell l 20051) . analyses of 
the evolution of the DF near a MBH have almost exclu- 
sively relied on the assumption that the mechanism through 
which stars exchange angular momentum and energy is dom- 
inated by uncorrelated two-bod y intera ctions (Chandrasekhar 
I 1943t see Binney & Tremaine 1987 for a more recent dis- 
cussion). This assumption is made in Fokker-Planck mod- 
els ( e.g. Ba hcall & WohfjUS [hereafter BW7 61; Bah call & 
WolfLL977| [hereafter BW77]; Cohn & Kulsrud 1978; Mur- 
phy, Cohn & Durisen Tl991l) . where the microscopic interac- 
tions are expressed by the diffusion coefficients, and in Monte 
Carlo simulations (e.g. Shapiro & Marchant 1979; Marchant 
& Shapiro l 19791.1 19801 Freitag & Benz l200lU 2002). Stars 
around MBHs are described as moving in the smooth aver- 
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age potential of the MBH and the stars, and the scattering by 
the fluctuating part of the potential is modeled as a hyper- 
bolic Keplerian interaction between a passing star and a test 
star. The scattering effects accumulate non-coherently in a 
random-walk fashion. 

The (non-resonant) relaxation time Tnr can be defined as 
the time Tg it takes for the negative specific energy £ = 
—v 2 /2 — ip of a typical star (hereafter "energy") to change 
by order unity. This is also the time Tj it takes for its spe- 
cific angular momentum J (hereafter "angular momentum") 
to change by an amount of order J c (£), the maximal (circu- 
lar orbit) angular momentum for that energy 3 . In a spherical 
potential J c = G[M, + N(> £)M*]/\/2£, where M. is the 
MBH mass, M* is the stellar mass and N(>£) is the number 
of stars with orbital energi es abov e £ (more bound than £). 
On Keplerian orbits J c = \fGM t a, where a is the semi-major 
axis. When relaxation is dominated by uncorrelated two-body 
interactions, the "non-resonant" relaxation time Tnr of stars 
of mass M* can be written in the Keplerian regime as 



Tmr = Ak — — 



M m 



P(a) 



M*J N(<a) 



(M.»M*), (1) 



where P — 2tt^J a 3 /(GM,) is the orbital period and Aa is 
a dimensionless constant which includes the Coulomb loga- 
rithm. We assume throughout a single mass stellar popula- 
tion; for numerical estimates we assume M* = 1 M Q . 

The assumption of uncorrelated two-body interactions is 
well-justified in many systems, such as globular clusters. 

3 Throughout this paper, "angular momentum relaxation time" means the 
time it takes until J is changed by order J c , rather than by order J. The time- 
scale for changes by order J < J c is shorter than the angular momentum 
relaxation time by a factor ( J/ J c ) 2 . 
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However, Rauch & Tremaine f 1996, hereafter RT96) showed 
that this does not hold for motion in potentials with certain 
symmetries, in particular for Keplerian motion in the poten- 
tial of a point mass where the orbits do not precess because 
of the 1 : 1 resonance between the radial and azimuthal fre- 
quencies. This is to a good approximation the case for stars 
orbiting close, but not too close to a MBH, since wider orbits 
precess due to the potential of the enclosed stellar mass, while 
tighter orbits precess due to General Relativity (GR). As long 
as the precession timescale t u (the time to change the argu- 
ment of the periapse u> by if) is much longer than the orbital 
period, the orbits effectively remain fixed in space over times 
P <C t -C t u and the interactions are correlated. The orbit- 
ing stars can then be represented by one-dimensional elliptical 
"wires" whose mass density varies along the wire in propor- 
tion to the time spent there (RT96). In this picture the interac- 
tion between two stars can be described as the torque between 
two wires. This results in mutual changes in both the direc- 
tion and the magnitude of the angular momenta. Following 
RT96, we denote such RR, which also changes the magnitude 
of J, scalar RR. The change grows coherently (oc t/t u ) on 

timescales t <C t w and non-coherently (oc yt/Tssd on longer 
timescales t ^> t u , where Trr is the "resonant relaxation" 
(RR) timescale for the angular momentum (Eq. [3}. Generally, 
Trr < T NR when M*N(< a) < M. (see Eq. Since the 
potential of the wires is stationary on time-scales <C t^, they 
do not exchange energy, and the energy relaxation timescale 
remains long, Tg ~ Tnr, even in the resonant regime. The 
mechanism of R R is rem iniscent of the Kozai mechanism in 
triple stars (Kozai 1962). 

A more restricted type of RR occurs in any spherical po- 
tential, where the vector J is conserved, but u> precesses. In 
that case the orbital rosette can be represented by a mass disk 
extending between the orbital periapse and apo-apse. The 
mutual torques exerted by such azimuthally symmetric disks 
change only the direction of J, but not its magnitude. This 
type of RR persists even in the presence of GR precession. 
Following RT96, we denote such RR, which changes only the 
direction of J, vector RR. The significance of the distinction 
between scalar and vector RR is that scalar RR can deflect 
stars into "loss-cone" orbits ( J < J c — * J ~ 0) where they 
fall into the MBH, either directly ("infall") or gradually ("in- 
spiral"), whereas vector RR can only randomize the orbital 
orientation, but cannot affect the loss rate. 

In this paper we explore the consequences of RR on the 
stellar DF and on the infall and inspiral rates. The mechanism 
of RR is briefly reviewed in ^2] The Galactic Center (GC) 
model is defined and the assumptions and approximations are 
discussed in ^3] In 50]We solve the Fokker-Planc k equ ation in 
energy and analyze the e ffects of RR on the DF ( 34.21 and on 
the stellar current ( ^4.3> . In ^5]we discuss the observational 
c onsequ ences of RR, such as the event rates of tida l disru ption 
( 35.1.H and gravitational wave (GW) emission f ^5. 1 .21 . We 
discuss our resul ts in the context of the young, non-relaxed 
stars in the GC ( 35.2> . We discuss and summarize our results 
infO 

2. RESONANT RELAXATION 

The resonant relaxation time Trr is estimated by evaluating 
A Ju, the coherent change in the magnitude of the specific 
angular momentum up to a time t u . The change A J u is then 
the step size ("mean free path") for the non-coherent growth 
of the angular momentum over times t > t u . Two nearby 
stars with semi-major axes a exert a mutual specific torque ~ 



GM+/a. Within a distance a from the MBH the net torque on 
a test star fluctuates away from zero as J ~ y/N(<a)GM+/ a 
and 

AJ W ~ jt w = y/N{<a){GMja)t„ . (2) 

For t > t u the torques on a particular star-wire become ran- 
dom, and the change in angular momentum grows in a random 
walk fashion with a timescale Trr ~ (J c / AJ u ) 2 t u , defined 
here as 

A / n r \ 2 



T RR = A 



RR 



N(>£) P\8) 
M 2 (>£) t u 



A 



RR 



M.Y P 2 {a) 



N(<a) \M*J i, 



(3) 

where /i = NM ir /(M,+ NM+), Arr is a numerical factor 
of order unity, to be determined by simulations, and the last 
approximate equality holds in the Keplerian regime. The TV- 
body simu lations of RT96 indicate that A RR = 3.56 for scalar 
RR (fJTTJ and A^ R = 0.31 for vector RR ( 353 (see their 
Eqs. 8, 15 and table 4 for /3 S and = (3~ 2 ). We adopt 

these values here. We verified these result, to within an order 
of unity, by carrying a few full small scale A-body simula- 
tions. Our simulations exhibited a substantial scatter in the 
value of ^4rr, as already noted by RT96. This introduces an 
uncertainty in quantitative estimates of the efficiency of RR. 
Further detailed analysis and simulations, outside the scope 
of this study, are needed to refine these numerical estimates. 

2.1. Scalar resonant relaxation 

Over most of the relevant phase space, the precession is due 
to the deviations from pure Keplerian motion caused by the 
potential of the extended stellar cluster ("mass precession"). 
This occurs on a timescale t u = tM, which assuming N(< 
a) AT <C M. and after averaging on J, can be expressed as 

M. 



t M = A M —, ^ P(a), (4) 

and where Am is a dimensionless constant. For simplicity 
we adopt here Am — 1, the value for circular orbits. The 
J-averaged RR timescale can then be written as 



M 
RR 



A^P(a) 



iRR 



A(<a)T r 



NR- 



(5) 



M. 

Since < T NR for small a where N(< a)M+ < A/., the RR 
rate of angular momentum relaxation is much higher than the 
rate of energy relaxation in the resonant regime. This qualita- 
tive analysis has been verified by detailed numerical A-body 
simulations by RT96 and by Rauch & Ingalls f 1998, hereafter 
RI98). 

For most of J-space, orbital precession is dominated by the 
mass of the stellar cluster and the RR relaxation timescale is 
well approximated by Trr ~ T^, which is a function of en- 
ergy only. However, when the orbital periapse is very close to 
the MBH, precession is dominated by GR effects ("GR pre- 
cession"). This is im portant for our analysis of the inspiral 
rate of GW sources ( 35.1.21 1. In this case the timescale for 
precession is given by t u = £qr> which is a strong function of 
J, 



*gr — — 



J 



3 V J T 



-'LSO 



P 



where 



J] 



AGM. 



LSO 



(6) 



(7) 
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is the angular momentum of the last stable orbit (LSO) for 
an orbit of specific kinetic energy <C c 2 . Generally, both 
precession mechanisms operate, and the scalar RR timescale 
t rr(£' J ) is S iven b y substituting l/t u = \l/t M ~ 1/*gr| 
in Eq. 0, where the opposite signs reflect the fact that mass 
precession is retrograde whereas GR precession is prograde. 
Thus, the scalar RR timescale is 



r, 



RR 



M. 



P\a) 



1 



1 



(8) 



N(<a) 

When tQR tjj. an d GR precession dominates, the RR 
timescale is (Eq.[3} 

Ig.^^Vf^V^. CD 



... P(a) 
I, J J iV(«z) 



The fact that scalar RR becomes much less efficient due 
to GR precession is crucial for the viability of GW emission 
from inspiraling stellar objects ("extreme mass-ratio inspiral 
sources", EMRIs), since it allows compact remnants to be 
very rapidly deflected to strongly relativistic orbits, but then 
stall "on the brink" and instead of falling directly into the 
MBH, inspiral into it gradu ally by the emission of GW. We 
return to this issue in 35.1.21 

We analyze below the dynamics of scalar RR and its effects 
on the stellar DF, both numerically ({J4j and by Monte-Carlo 
simulations These two approaches require different for- 
mulations of the RR timescale. The numerical analysis is car- 
ried out in the ID 5-space, where the J-dependences are ab- 
sorbed in the J-averaged terms. For this purpose we use the 
relation d( J 2 )/ J c 2 = di/T RR (£, J) (Eqs. [S> to define the J- 
averaged time it takes a star to random- walk from J = J c (£) 
to the loss-cone J = J; c as 



rr(£) - ~p 



dJ 2 T£ R (£,J). 



(10) 



Since at low energies T RR oc P, it initially decreases 
with £ toward the MBH. However, at very large energies 
near the MBH, where GR precession becomes important 
(JhSO / Jc(£) —> 1), 2r R increases again (Fig. Q. The Monte 
Carlo simulations are carried out in 2D (£,J)-space. For this 
purpose we define a general expression for the scalar angular 
momentum relaxation time, Tf, which applies in any regime 
(Eqs. [T]|3l|8l see also footnote^). 



T}(£,J) 



1 



1 



T NR (£) Trr(£, J) 



(11) 



2.2. Vector resonant relaxation 



Vector RR grows coherently (oc t) on timescales t -C t v , 
where t„ is the timescale for a change of order unity in the 
total gravitational potential ip caused by the changes in the 
stellar potential due to the realignment of the stars as they 
rotate by n on their orbits, 



t -2--A 



N X ' 2 P 



2 



Mm P 



(12) 



where A v is a dimensionless constant of order unity, and 
where the last approximate equality holds for NM+ <C M,. 
For simplicity we adopt here A v = 1. In analogy to scalar RR 
(Eq. |3, the maximal coherent change in J is 



I AX, 



Jt u 



Jc 



(13) 



that is, J rotates by an angle 0(1) already during the coherent 
phase. On timescales t ^S> t v , |AJ V | cannot grow larger, as 
it already reached its maximal possible value, but the orbital 
inclination angle is continuously randomized non-coherently 
(oc t 1 / 2 ) on the vector RR timescale (Eq. [3}, 



1 RR 



2A 



v N ^ 2 (>£) 



RR' 



P{£) ~ 2A 



RR 



Mm 

m7 



P(a) 



W 1 /2(< a )' 

(14) 

where the last approximate equality holds for iVM* <c M m ■ 

Note that while the torques driving scalar and vector RR 
are the same, vector RR is much more efficient than scalar 
RR, T^ R <C ?r R , due to the much longer coherence time 

ttp ~ N l ' 2 tM ^S> tj\f- Furthermore, vector RR proceeds ir- 
respective of any precession mechanisms that limit the effi- 
ciency of scalar RR. 

3. A STELLAR CLUSTER WITH A MASSIVE BLACK HOLE: MODEL 
AND ASSUMPTIONS 

In this section we state our assumptions, approximations 
and definitions, many of which are commonly used in the 
analysis of stellar systems with MBHs. Since the GC is the 
best studied case of a stellar cluster with a MBH, and since 
it is representative of typical LISA targets, we scale the pa- 
rameters to values appropriate for this system. Our goal is to 
find the stellar DF in presence of RR. A full solution in (£ , J)- 
space is complex and is not attempted here. For simplicity, the 
analysis carried out here is in £-space only and it assumes an 
underlying spherical symmetry. The J-dependent loss-cone 
terms are approximated by their J-averaged effective terms, 
which are functions of £ only. 

3.1. Distribution in energy space 

Recent observations show that a tight empirical relation ex- 
ists between the mass of MBHs, M # , and the velocity disper- 
sion a of their host bu lge or galaxy (Ferrarese & Merriti 2000; 
Gebhardt et al J"200Ch . 

Mm = 1.3 x 1O*M (^pr) 4 , (15) 

where a is me asured at the galaxy's effective radius 
(Tremaine et al. 2002). Both the exponent and the pre- 
factor have small uncertainties, which we ignore here. The 
Galaxy obeys the M,-a relation, with AT. ~3x 10 6 Af Q , and 
a ~ 75 km s ~ 1 (Tremaine et al J 20021) . 

The MBH dominates the dynamics of stars within its 
"Bondi radius", or radius of influence, 



CM, 

r h = — — = 2 pc 



Mm 



= 2pc 



Mm 



3 x 10 6 M Q 

1/2 



3 X W>Mr. 



(r5k^)~ 2(16) 

(17) 



where in the last step we used the M,-a relation. Within r;, 
the system is approximately Keplerian, and deviations from 
Keplerian motion become important only for times t ^> P, 
where P(£) = 2tt x /(GM.) 2 /8£ 3 is the orbital period and 
the energy is £ = GAl m /r — v 2 /2. 

We assume a single mass population with an underlying 
spherically symmetric DF, so that /(x, v) = f(£) (stars per 
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phase space volume d 3 xd 3 v). Following BW76 we assume 
that for r > rh the DF is Maxwellian, 

m = j2^ e£ ' al {r>rhh (18) 

where rift w 4 x 10 4 pc~ 3 is the number density at ru in the 
GC assuming a mean mass of 1 M Q (Genzel et al. 20Q3J) and 
where 07,. = 75 km s _1 is the ID velocity dispersion at r^. For 
Keplerian orbits, the density n(£) (stars per d£) is 

n(£) = ir 3 V2(GM.) 3 £- 5 / 2 f(£) , (19) 
the density n(r) (stars per d 3 x) is 

r GM./r fnT7 

n(r)=AV2n d£f(£)J '--£, (20) 

J — 00 » ' 

and the density n(a) (stars per da) is 

n(a)=4, 3 (GM.) 3 /V/ 2 /(^) . (21) 

Two stellar components contribute to the local stellar den- 
sity at radius r. One is the tightly bound stars with energy 
£ ~ 0(GM,/r), or a ~ 0(7"). The other is the unbound, or 
marginally bound stars with energy £ <C GM,/r (or a^> r) 
and high eccentricities, which spend only a small fraction of 
their orbit inside r. When unbound stars dominate the local 
population, for example because RR or stellar collisions have 
destroyed most of the tightly bound stars there, the density 
profile is given by 



n u (r)=4V27r J d£f(£)J^- - £ cx r" 1 / 2 . (22) 

We express the local relaxation time, which generally de- 
pends on the nature of the relaxation process and on radius, 
in terms of a reference time Th, which we define as the NR 
relaxation time at 77,, 

h 16 V n n h (GM*) 2 In A 

( M \ 3/4 

= 6xloV M^) ■ < 23 > 

where the last equality assumed the M»-a relation and the 
stellar density in the GC. The Coulomb factor is A = 
fmaxAmm = M./M* (BW76), where r max ~ r is the max- 
imal impact parameter for perturbations by stars interior to r, 
and r m ; n is the minimal impact parameter for small deflec- 
tions, r min = GM±/v 2 = (M*/M.)r. The local NR relax- 
ation time may be quite different from Th when the density 
distribution of the system strongly deviates from the steady 
state configuration n(r) cx r~ 7 / 4 ( 34.2> . 

3.2. The loss-cone 

At high enough energy the DF vanishes, /(£ > £d) = 0, 
because stellar objects cannot survive close to the MBH. The 
value of £0 depends on the process that destroys the stars. If 
they are compact objects, they fall directly into the MBH or 
inspiral into it by GW emission. If they are main sequence 
stars, they are destroyed by disruptive stellar collisions in the 



high density cusp around the MBH, or are tidally captured 
and heated until disruption, or are tidally disrupted when their 
orbital periapse r p falls below the tidal disruption radius, r t ~ 

(M./M*) 1/3 R±. 

Stars are destroyed by the MBH either by scattering or de- 
caying in £-space to the point where £ > £jj, or by scat- 
tering in J-space to the point where J < J; c = \f2GM,r t 
(or r p < r t ). Scattering in J-space is by f ar mor e efficient 
(Frank & Rees 1976; Lightman and Shapiro 19771 see foot- 
note|3jl, and the tidal disruption rate is dominated by stars on 
low-£ and low- J orbits. Stars enter the loss-cone, defined by 
J < Jic, in two regimes of phase-space, depending on their 
orbital energy, which determines the ratio between the angu- 
lar momentum change per orbit, A J(£), and J; c . For small 
energies (long periods), A J ^ J/ c . In this "full loss-cone" 
(or "kick") regime, the DF remains essentially unmodified by 
the existence of a loss-cone. For high energies (short periods) 
near the MBH, A J<C Ji c - In this "empty loss-cone" (or "dif- 
fusive") regime stars slowly diffuse into the loss-cone and are 
then promptly destroyed on a dynamical time, and the DF is 
modified by the existence of the loss-cone. The stellar density 
vanishes inside the loss-cone, n(£ , J < J; c ) = 0, and gradu- 
ally falls to zero to ward it, n(J) cx Jln(J/J/ c ) for J > Ji c 
(Lightman & Shapir o! 9771) . This also modifies somewhat the 
DF in £-space. Note that scalar RR is typically efficient only 
deep in the empty loss-cone regime, where mass precession is 
negligible. 

There is a critical energy where AJ(£) = J; c , that de- 
marcates the transition between the empty and full loss-cone 
regimes. The size of the loss-cone, and therefore the critical 
energy, depends on the nature of the loss-process. The critical 
energy for prompt stellar disruption £ v , corresponds to a dis- 
tance scale of ~ r/j where £ = £ h (Eq. I I6I 1. Most of the contri- 
bution to the stellar destruction rate is from stars with £>£ p . 
Similarly, there is a critical energy for stellar destruction by 
GW inspiral, £gw- GW inspiral takes much longer than direct 
disruption. Unless the star starts out on a short-period (high- 
£) orbit, it will be scattered again directly into the MBH, or 
to a high- J orbit where GW dissipation is negligible. Thus, 
£cw^£p (Alexander & Hopman 200jt. 

We adopt below for numerical calculations the values of the 
critical energies in the GC, where for consistency with our 
assumed single mass population, we consider only prompt 
disruption and GW inspiral at the last stable orbit, J/ c = 
4GM,/c. In the GC, the semi-major axes corresponding to 
these critical energies are a p = 0.27 pc and acw = 0.02 pc 
(see Fig. |6j. Thus, our model describes a population of com- 
pact objects, and our predicted prompt disruption rates will 
be slightly lower than the actual tidal disruption rates, which 
occur at the somewhat larger J; c = y/2GM,rt and thus have 
a lower critical energy, closer to £ h. 

A complete treatment of the loss-cone problem involves the 
solution of the Fokker-Planck equation in (£, J) -space (Cohn 
& Kulsrud 1978). However, an approximate solution can be 
obtained by solving the Fokker-Planck equation in £-space 
only, while accounting for the loss-cone by adding sink terms 
to the equation (BW77). This is the approach we adopt here. 
The tidal disruption rate due to NR relaxation has been dis- 
cussed extensively in the literature (see e.g. Lightman and 
Shapiro Tl977t Frank & Rees Tl976l Cohn & KulsrudQ97! 
Syer & Ulmer 1999; Magorrian & T remain e l 1999t Alexan- 
der^ Hopma j2003i Merritt & Pood 2004 Wang & Merrit t 
I 20041 Hopman & Alexander 20051 Baumgardt et al. 1 20051) . 
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The differential NR loss-cone diffusion rate (stars per d£dt) 
is estimated to be of the order 



i NR (£) ~ n(£)/T NR (£) 



(24) 



neglecting here the weak logarithmic J-dependence that ex- 
presses the depletion of phase space near the loss-cone. When 
RR is the dominant relaxation process, then for timescales 
3> t u the differential loss-cone refilling rate is 4 



T RR (f ) ~ n{£ )/fi R (£ ) 



(25) 



Note that vector RR does not enter into the Fokker Planck 
equation and does not play a role in the diffusion in £-space 
since the system is assumed to be spherically symmetric, and 
vector RR does not change J. 

4. DENSITY PROFILE NEAR A MASSIVE BLACK HOLE 
4.1. Fokker-Planck energy equation with resonant relaxation 
Following BW76, we define the dimensionless quantities 

.9 = (2 7 ra2) 3/2 n / ; 1 /; x = £/a 2 h : r = t/T h , (26) 
and write the Fokker-Planck equation 5 as 



dg{x,r) = 5/2 3Q(x,t) 
dr dx 



Rnr{x,t) - xRrr{x,t), 
(27) 



where the current Q(x), defined to be the net rate at which 
stars flow to energies > x, is given by (BW76) 



Q(x,t) 



dy 



dy 

x {max(x, y)} 3 ^ 2 . 



dg{x,- 
dx 



(28) 



The sink terms -Rnr and Rrr represent the J-averaged dif- 
ferential loss rate of stars per energy interval, by NR and RR 
respectively, due the existence of a loss-cone in J-space. This 
is an approximate substitute for the full 2 + 1 treatment. 

The NR sink term R^ R (x, t) is expressed by the form (Eq. 
IB 1 l in appendix B) proposed by BW77, which smoothly inter- 
polates between the empty and the full loss-cone regimes and 
approximates the logarithmic depletion of phase space near 
the loss-cone boundary ( §3.2> . In the empty loss-cone regime, 
which is of interest here, 



RNR(x,T)(xg(x,T) 2 /log(x D /4x) 



(29) 



This quadratic behavior arises since a fraction g(x, r) of the 
stars is accreted per local relaxation time, which is itself pro- 
portional to 1 /g(x, t). 
The RR sink term Rrr is expressed by (Eas.ElEBl 



Rrr(x, t) = g(x, t) /trr(x) 



(30) 



where trr(x) = Tp R (£)/J\ is the dimensionless J-averaged 
RR time, defined in Eq. fllOt . 

4 Eq. 1251 differs from the corresponding term in RT96, which appli es on ly 
for timescales <g tw, and an isotropic DF. The rate predicted by Eq. 1251 is 
in good agreement with the values obtained by the iV-body simulations of 
RI98. 

5 Equation 1271 is written in the form of a particle conservation relation, 
which, as shown by BW76, is equivalent to the usual Fokker-Planck form 
df/dt = —D x df/dx + D xx d f / dx 2 , where D x and D xx are the diffu- 
sion coefficients. 



The RR efficiency factor \ m Eq. (I27> parametrizes the 
uncertainties in the efficiency of RR and in the approxima- 
tions involved in Eq. (I30i . for example the effect of the par- 
tial depletion of phase-space outside the loss-cone (note that 
the depletion is small, see RI98 fig. [2a]). For simplicity, we 
assume that \ is a constant. We expect and show 

below that the choice x = 1 yields results that are consistent 
with those of RT96 and RI98, but we also explore RR for a 
range of \ values to determine the robustness of our results. 
The RR sink term becomes negligible compared to the NR 
sink term for lower energies, were Trr ^> Th (Fig. Q. 

We follow BW76 by assuming as boundary conditions 



g(x<0,r) = e* 



g(x = x D ,r) = 0. 



(31) 



The first boundary condition for x < expresses the assump- 
tion of a thermal reservoir at the radius where the stars are no 
longer bound to the MBH (Eq. II 81 . In reality, they are bound 
to the total enclosed mass of the MBH and stars, however the 
BW76 treatment neglects the stellar mass and assumes Ke- 
plerian motion around the MBH. The second expresses the 
existence of a mass-sink at energy xd ■ Because of the large 
dynamical range in the problem, it is more convenient for nu- 
merical purposes to express Eqs. ( 12713 1> in logarithmic en- 
ergy (see appendixlAV 

4.2. The effects ofRR on the stellar distribution function 

Here we explore the effect of RR on the steady state DF. We 
begin by checking the numerical convergence of Eq. Mil to 
the BW76 solution, f(£)<x£ p withp= 1/4 for £ h -^£<^£ D , 
which is obtained under the assumption that the £ -dependent 
J-averaged sink terms can be neglected, i?NR = ^rr = 0, and 
that the disruption of stars by the MBH can be expressed by 
the boundary condition f(£ > £d) = 0. Figure Q confirms 
that the DF indeed converges to this solution on a relaxation 
timescale, and we also verified that this holds irrespective of 
the initial DF assumed. Rapid conversion was also confirmed 
by the iV-body simulations of Merritt & Szell If 20051) . The 
BW76 DF corresponds to a power-law density profile n(r) cx 
r~ a with a =p + 3/2 = 7/4 (eauationl20l Since T NR (£) cx 
P{£)/N(> £), and T NR (£ = a\) = T h , the implied local 
NR time is 



7nr(£) = T h (£/a 2 h ) 



3/2- 



(32) 



We now add back the NR and RR sink terms and solve the 
full problem. Figure (0 shows the results of the steady state 
solution of equation (127 \ for several values of \. The \ = 
case (NR only) reproduces the solution studied by BW77, 
who showed that the presence of the NR sink term does not 
change the power-law behavior and index of the DF. We find, 
as expected, that for low energies where RR is negligible, nei- 
ther does the presence of the RR sink term. At higher ener- 
gies, the response of the DF depends on the efficiency of RR. 
For x ^ 10, RR is so efficient that it exponentially depletes 
the high energy end of the DF. However, for the likely effi- 
ciency factor x = 1) GR precession at high energies limits the 
efficiency of RR (cf Fig. 0, and the DF is not completely de- 
pleted, but continues to rise toward £ — > £r> after a moderate 
initial drop. 

4.3. 77ie effects ofRR on the stellar current 

The solution of the Fokker-Planck equation determines the 
stellar current Q(x) into the MBH. This in turn determines 
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FIG. 1. — The evolution of the stellar DF to a steady-state f 1 / 4 cusp 
(BW76) by ^-diffusion only (no loss-cone effects: i?NR = Rrr = 0). The 
dimensionless sink energy is at xo = 10 5 , approximately appropriate for 
the tidal radius in the GC. Lines show g(x, r) at intervals of At = 0.07, 
with the lowest curve indicating the initial DF (arbitrarily assumed here to be 
g(x, 0) oc x 0-1 ). The system was integrated up to 21),; the steady state DF 
is reached after ~ 0.8T^. For comparison, an x 1 / 4 power-law is also shown. 
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FIG. 2. — Growth of a cusp in presence of RR, with the same initial 
conditions as in Fig. Q in the presence of the NR and RR sink terms for RR 
efficiency factors of \ = 0, 1, 3, 10. For RR efficiency factor \ ^3, RR leads 
to an exponential cutoff of the DF at high energies. 



the rates and modes of stellar capture by the MBH. Since 
most of the contribution to prompt infall events comes from 
o rbits with energy x>x p , and to GW inspiral from x > xqw 
( SI3.2I ^5}, the current at these critical energies determines the 
prompt infall and gravitat iona l inspiral event rates. In the GC, 
x p = 2.7 and xqw = 58 ( 33. 2> . 
The rate at which stars flow toward the MBH is given by 

I(£,t) = I Q(x,t), (33) 
where the dimensional current scale is set by 

_ Svr 2 3 (GMJ 2 In An h ( N h \ 

/o = uf hnh — 3 ° UJ • (34) 

For non-equilibrium DFs, Q ~ 1, and there is a large cur- 
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FIG. 3. — The steady state current Q as function of dimensionless energy x, 
for various choices of the RR efficiency parameter \. Solid line: \ = (NR 
only). Dashed line: x = !• Dot-dashed line: \ = 10- The capture rates are 
determined by the value o f the current a t the critical energies (shown here at 
values typical for the GC, 315.1.11 85.1.21 . The dotted line is the ratio Q\j Qo 
between the currents of the \ = 1 and \ = cases. At iqw. Qi/Qo ~ 8, 
so the GW inspiral rate is dominated by RR. At x p , Qi/Qo ~ 3 so prompt 
infalls are somewhat enhanced by RR, as found by RT96 and RI98. When 
X = 10, RR is so efficient that it evacuates a large fraction of the stars at 
energies higher than ig^. and me GW inspiral rate is suppressed relative to 
the x = l, whereas the prompt infall rate is further enhanced. 



rent I ~ Iq. The "sink-less" BW76 solution ( §4.2t has a 
steady state "zero-current" (0 < Q -C 1) solution that is £- 
independent and is strongly suppressed by the bottle-neck at 
£d- The presence of the loss-cone, as expressed by the NR 
and RR sink terms, modifies the DF at high energies and shifts 
the effective high-energy boundary to lower energies. As a 
result, Q increases well above the zero-current solution and 
becomes ^-dependent. 

Figure l|3} shows the current Q(x) for several values of the 
RR efficiency factor x- The accelerated relaxation due to RR 
increases the current above the values induced by NR only. 
For x — 1 me largest enhancement is attained at about the 
critical energy for GW. When x ^ 10, RR is so efficient that 
the current depletes the DF already at energies x < xqw ■ As a 
consequence, the current drops to zero at x > xgw, and there- 
fore so does the GW inspiral event rate (which is oc Q(o;gw)> 
see eq. 1371 ). Figure @ shows separately the contributions 
of the NR and RR sink terms in Eq. J27I to the total capture 
rate, and demonstrates that far from the MBH (low-£), where 
most tidally disrupted stars originate, NR dominates the cap- 
ture rate, whereas at closer distances (high-£), where GW in- 
spiral stars originate, RR dominates the capture rate (see also 
Fig-Qbelow). 

5. OBSERVABLE IMPLICATIONS OF RESONANT RELAXATION 

Most previous estimates of the loss rates did not take into 
account RR. Here we show that RR can strongly influence the 
GW inspiral rate, and may play an important role in determin- 
ing the dynamical structures observed in the inner GC. 

5.1. Infall and inspiral processes 

5.1.1. Prompt infall 

Figure Q shows the change in the estimated rates of 
prompt infall and inspiral due to the inclusion of RR, as func- 
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FIG. 4. — The differential capture rates x 5 / 2 Rnh (solid line) and 
a; 5 / 2 -RriR with \ = 1 (dashed line) as function of dimensionless energy 
x. At high energies near iqw. where GW sources originate, captures are 
dominated by RR. Prompt infall is dominated by NR in the region near x p , 
in the transition between the full and empty loss-cone regimes. 
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FIG. 5. — The x-dependence of the enhancement of the prompt infall and 
GW inspiral event rates by RR over NR, as expressed by the stellar current at 
the critical energy x c , Q(x c ; x)/Q{x c \ 0), for x c = x p , xqw- 



tion of the RR efficiency parameter \- RR does not signifi- 
cantly increases the rate of prompt infall events. This result 
has already been demonstrated by the A-body simulations of 
RT96 and RI98, and is explained by the fact that most tidally 
disrupted stars originate from low-£ orbits wi th £ ~ 0(£h) 
(Lightman & Shapiro Tl977t Baumgardt et al. I 20 05). where 
mass precession suppresses scalar RR. Our analysis is com- 
plementary to that of RT96 and RJ98 in that they took into 
account relaxation in J-space, but not in 5-space, whereas 
ours explicitly solves self-consistently the DF and star current 
in £-space, but does not deal with the flow in J-space. Our re- 
sults confirm the conclusion that for the likely value of the RR 
efficiency factor x = 1, the tidal disruption rate is enhanced by 
RR by only a factor of ~ 3 (Figs[3][5j. We note, however, that 
for higher values of x> the tidal disruption rate may increase 
significantly. 



5.1.2. Gravitational wave inspiral 

The detection of GW emission from compact remnants in- 
spiraling into MBHs (EMRIs), is one of the major goals of 
the planned space -borne LISA GW detector 6 . The predicted 
event rate per galaxy is very uncertain, as estimates vary in 
the ra nge Gy r -1 < Tgw My r -1 per galax y (e.g. Hils & 
Bended 19951 Sigurdsson & Ree s! 19971 Iv ano M 2002b Freitag 
2001 2003; Hopman & Alexande ?2005l) . 

Only stars originating from tightly bound orbits with ener- 
gies £ > £gw can complete their inspiral and reach a high- 
frequency orbit observable by LISA, without being scattered 
prematurely into the MBH or to a wider orbit. An approxi- 
mate upper-limit for the semi-major axis from which inspiral 
is possible can be obtained by equating the inspiral time tew 
with (1 — e)TNR, the N R time scale for scattering by A J ~ J, 
(Hopman & Alexande J^OOH) . 

flGW H ^ ) ' El = 3x213 M. ' (35) 

(the corresponding critical energy is £gw = GM,/2acw)- 
There are three phases in the orbital evolution of a star that 
ends up as a GW source. The first is a scattering-dominated 
phase in J-space, where energy losses are negligible. This 
is followed by a transition phase, when J ~ few x Jlso, 
GW dissipation becomes significant, and the inspiral time 
is comparable to the scattering time. Finally, the star spi- 
rals in by GW emission so rapidly that it effectively decou- 
ples from the gravitational perturbations of the background 
stars. The value of acw is determined in the second, transition 
phase by the interplay between angular momentum relaxation 
and GW dissipation. The J-values typical of the transition 
phase are also those where GR-precession quenches RR, at 
J < Jq ~ fewxJLso (Eq. |6j. As a result, ogw is still deter- 
mined by NR through Eq. 1351 even when RR dominates the 
dynamics at J^> Jlso- 

To verify that RR does not significantly affect the transition 
phase of inspiral and the value of the critical energy, we fol- 
lowed the approac h used by Hils & Bender f 1995) and Hop- 
man & Alexander (2005), who performed Monte Carlo (MC) 
simulations of the scattering of stars in J-space. We compare 
the case w here Tj = Tnr (as was ass umed in Hils & Ben- 
der 1995; Hopman & Alexander 2005) to the case where Tj 
includes RR (Eq. 111b . Figure shows our results in terms 
of the inspiral fraction S(a), the fraction of MC experiments 
that end with the star spiraling in, relative to the total, which 
also includes the experiments that end with the star plung- 
ing directly into the MBH without emitting a detectable GW 
signal (the third possibility of diffusion to lower energies is 
not relevant here since only J-scattering is considered). The 
critical semi-major axis corresponds to S(clgw) ~ 0.5. Fig- 
ure shows that RR has little effect on the value of acw, 
which is completely determined by NR. Furthermore, the out- 
come is not sensitive to the exact value of Jq where RR is 
quenched, as is demonstrated by the fact that oqw is only 
slightly decreased when the strength of GR precession is ar- 
tificially lowered by a factor of 10 (by increasing x t0 10), 
which is equivalent to assuming a smaller value of Jq. 

Because S(a) is well approximated by a step function at 
bgw (or equivalently, S{£) is a step function at £gwX the 
inspiral event rate per galaxy without RR can be estimated by 

6 http://lisa.jpl.nasa.gov/ 
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FIG. 6. — The fraction of stars S(a) that produce a detectable GW signal 
out of all stars consumed by the MBH, as function of the initial semi-major 
axis. The solid line is for NR only; the dots for NR and RR (including GR 
quenching near J; c ). The two plots are almost identical, proving that RR 
does not affect the inspiral probability. The dashed line is for NR and RR, 
but with an artificially decreased GR precession (by a factor 10). In this case 
S(a) drops at slightly smaller values of a. The small difference implies that 
S(a) is not sensitive to the exact value of Jq where RR is quenched. 
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FIG. 7. — Stellar components, timescales and distance scales in the GC 
(equating a ~ r). The NR timescale Tnr (top straight line); the averaged 
scalar RR timescale T^ R , estimated for 1 Mq stars (top curved line) and 
10 Mq stars (bottom curved line); the vector RR timescale Tp R (bottom 
straight line); the position and estimated age of the young stellar rings in the 
GC (filled rectangle in the bottom right); the position and age of the S-stars 
if they were born with the disks (empty rectangle in the bottom left); the 
position and maximal lifespan of the S-stars (filled rectangle in the middle 
left); the position and age of the dynamically relaxed old red giants (filled 
rectangle in the top right); and the position and age of the reservoir of GW 
inspiral sources, where the age is the progenitor's age or the time to sink to 
the center (empty rectangle in the top left). 



(Hopman & Alexandei 2005) 



r 



GW 



d£ n(£) 



few MJc/Jlso)T nr (£) 



(36) 



where f s is the population fraction of the specific GW sources 
under consideration (e.g. white dwarfs). 

Figure shows that the GW inspiral event rate reaches its 
maximal enhancement, of about an order of magnitude, for 
X — 1 and is enhanced over a wide range of values. Only 
when x ^ 10 is the RR depletion of high energy GW inspiral 
candidates with x > xqw so strong (Fig. |2} that the GW in- 
spiral rate falls below that predicted for NR only. This deple- 
tion can be expressed by introducing an effective high-energy 
cutoff £rr <C £d, such that n{£) ~ for £ > £r,r. Two 
very different situations can occur, depending on whether the 
critical energy £qw is larger or smaller than the cutoff energy 
£rr- 

If £gw ^S> £rr, the DF vanishes for all energies where stars 
could complete their inspiral. In this case the GW event rate 
is vanishingly small. Fortunately for the prospects of EMRI 
detection, it appears likely that this is not the case (see Figs. 
E1EH3- When £ G w<£r, the GW event rate is given by 



d£n{£) 



= fsIoQ(xGw) , (37) 



'few ^Rr(£) 

where the last equality follows from Eq. (I27> in steady state. 

The Galactic MBH is the prototypical LISA target (the or- 
bital frequencies around more massive MBHs are below the 
LISA sensitivity band). Applying our result to the GC, we 
estimate that RR increases the predicted LISA event rate per 
galaxy by up to a factor (Fig. [3} 



rf Y = Qrr(zgw) 
Tnr <3nr(£gw) 



(38) 



over previous predictions that assumed NR only. One caveat 
is that our analysis neglects the effects of mass segregation 



in a multi-mass population. For example, low-mass remnants 
will be driven by mass segregation to a flatter density distri- 
bution, which will decrease their contribution to the inspiral 
event rate (Hopman & Alexander 20051 1 20061) . and the op- 
posite will happen for high mass objects. Here we do not 
consider a multi-mass population. 

5.2. Resonant relaxation in the Galactic center 

At a distance of — 8 kpc (Reid & Bruthaler QOQ^) . the 
Galactic MBH is the closest and most accessible MBH. As- 
trometric and radial velocity measurements of stars closely or- 
biting the MBH indica te that its mas s is M. ~ (3-4) x 10 6 M© 
(Scho del et al. I 20021 Ghez et al. I 20031 Eisenhauer et al. 
2005; for a review see Alexander 2005). Observations of 
the GC provide the most detailed information available about 
stars close to a MBH, where the orbits are Keplerian. The 
Galactic MBH offers, in principle, the best chances of con- 
fronting RR with observations. We consider here the impli- 
cations of RR for several observed and predicted features of 
the GC: the young coherent star disks, the central randomized 
cluster and the relaxed old giants; the orbital parameters of the 
innermost stars; the stellar density distribution; and the hy- 
pothesized dense stellar cluster of stellar black holes (SBHs). 

Figure |7) compares the distance scales and the ages or 
lifespans of the various dynamical structures and components 
in the inner pc of the GC with the relaxation timescales. The 
NR timescale in the GC, which is roughly independent of ra- 
dius, is Tnr ~ fewx 10 9 yr (Eq. [Q. The scalar RR relax- 
ation time Tft R is shown for x — 1, 10, or, equivalently, for 
M* = 1, 10 M Q with x — 1- At large radii the RR time de- 
creases towards the center, but for small radii, where GR pre- 
cession becomes significant, it increases again (see 50. The 
vector RR timescale T^ R , in contrast, decreases unquenched 
with decreasing radius. Structures whose estimated age ex- 
ceeds these relaxation timescales must be relaxed. Structures 
whose lifespan exceeds the relaxation timescales may be af- 
fected, unless we are observing them at an atypical time soon 
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after they were created. As is further discussed below, RR can 
naturally explain some of the systematic differences between 
the various dynamical components in the GC. 

The observed deviations from spherical symmetry of the 
potential in the GC do not affect RR. As an example, one may 
consider a ring of mass M r and radius r. Such a perturbation 
to the potential would change orbits with semi-major axis a on 
the Kozai time fKoza H^I) t K ~ 2ir{M./M r ){r/a) 3 P{a), 
which is larger than tj^ by a factor of order / tj^ ~ 2ttN(< 
a)M^/M r (r/af < 2wN(< r)M*/M r . This implies that 
only if the mass of the ring is comparable to the enclosed 
mass, such a perturbation to sphericity will be important. This 
is not the case in the GC. For example, the heaviest young stel- 
lar disk has a mass M r < 10 4 M Q (Nayakshin et al.[2006), 
giving tx/tM > 400 for r — 0.4 pc. 

5.2.1. RR and the young stars in the Galactic center 

Two distinct young stellar populations in the GC may be of 
particular relevance for testing RR. At distances of 0.04-0.5 
pc from the MBH there are about ~ 70 young massive OB 
stars (M* > 10 M , lifespan of = 6 ± 2 Myr), which are 
distributed in two nearly pe rpendic ular, tangentially rotating 
disks (Lev in & B elobodorov 2003; Genzel et al. 20031 Pau- 
mard et al. l 20051) . It appears that these stars were formed by 
the fragmentation of gas disks (Levin & Belobodorov 2003; 
Levi nl 20031 Nayakshin & Cuadra l 2004l Nayakshin & Sun- 
yaev[_2005; Nayakshin 2006). This young population co- 
exists with a relaxed population of long-lived evolved giants 
(U > 100 Myr, M* < 8M Q ; Genzel, Hollenbach & Townes 
1994). Inside the inner 0.04 pc the population changes. There 
is no evidence for old stars, and the young stars there (the 
"S-stars") are main-sequence B-stars (M + < 1 5 Mq, lifes- 
pans ofl0 7 < t* <2 x 10 8 yr; Ghez et al. l 20031 Eisenhauer 
et al.Q005) on randomly oriented orbits with a random (ther- 
mal) J-distribution. The orbital solutions obtained for a few 
of the S-stars indicate that they are tightly bound to the MBH 
(Ghez et al. l 20031 Eisenhauer et al. l 20051) . There is to date 
no satisfactory explanation for the presence of the S-stars so 
close to the MBH (see Alexander 2005| for a review). One 
possibility is that they originated far from the MBH and were 
captured near it by dynamical exch ange in teractions (Gould 
& Quille d 20031 Alexander & Livi d 20041) . In that case only 
their lifespan can be confidently determined, but not their ac- 
tual age. Another possibility is that they also originated in the 
star disks. In that case their age is the same as that of the disks 
and is much shorter than their lifespan. 

The existence of coherent dynamical structures in the GC 
constrains the relaxation processes on these distance scales, 
since the relaxation timescales must be longer than the struc- 
ture age i* to avoid randomizing it. Figure |7} shows that the 
observed systematic trends in the spatial distribution, age and 
state of relaxation of the different stellar components of the 
GC are consistent, and perhaps even caused by RR. The star 
disks are young enough to retain their structure up to their 
inner edge at 0.04 pc, where t* ~ T RR and vector RR can 
randomize the disk. It is tempting to explain the S-stars as 
originally being the inner part of the disk. However, RR alone 
cannot explain why the S-stars are systematically less massive 
than the disk stars. Regardless of their origin, their random or- 
bits are consistent with the effect of RR. Vector RR can also 
explain why the evolved red giants beyond 0.04 pc 2 in par- 
ticular the more massive ones with min(X , NR, ^j^j^) are 
relaxed, since T RR < t+ out to ~ 1 pc. Note that scalar RR is in 
itself not efficient enough to drain the S-stars, but it can play 



a role in their non-zero eccentricities (see also Levin 2006). 

Rapid RR changes the orbital eccentricities. On the 
timescales relevant for observations (t <C t u , Eq. 0}, the 
change grows linearly and could be detected, in princi- 
ple, by precision astrometry. However, the effect is small. 
In the linear regime the shortest time-scale for changes 
in J is determined by scalar RR. The fractional change 
per orbit is AJ P /J C ~ y/N(< a)(M*/M.) (Eq. EJi, 
which corresponds to a change in eccentricity of Aep ~ 
y/N(<a)(M i ,/M,)y/l — e 2 /e. For the best measured or- 
bit (th e star S2 with e = 0.8760 ± 0.0072. Eisenhauer et al. 
2005) the predicted change per orbit is Aep ~ 10~ 4 (assum- 
ing M* = 10 M , a = 0.01 pc and AT(< a) = 2500). This is 
beyond the current astrometric precision. 

5.2.2. RR and the stellar distribution in the Galactic center 

Figure (|8j shows the density distribution n(r) (Eas. 12011221 
and the distribution of semi-major axis a~ 2 n(a) (Eq. 12 11 1 for 
a model of the GC with different values of x- The central 
depletion in the distribution of the semi-major axis (or equiv- 
alently, energy) does not appear in the space density profile, 
because unbound and weakly bound stars on eccentric orbits 
spend a fraction of their orbit in the center. The effect of RR 
on the density distribution only appears as a gradual flatten- 
ing of the cusp slope to n u {r) cx r -1 / 2 (Eq. I22> . The effect 
on the observed projected surface density distribution will be 
even less noticeable. RR is predicted to flatten the density dis- 
tribution only at < few x 10 -3 pc, for x = 1- Farther out the 
distribution is virtually identical with the NR solution. Since 
there are no observational data on the density distribution for 
distances <5x 10~ 3 pc (Schodel et al. 2005), it is presently 
not possible to test this aspect of RR in the GC. However, if 
X 3> 1, the space density will flatten farther out, on a length- 
scale that may be observationally accessible. It should be em- 
phasized that any detailed comparisons between the observed 
stellar distribution and models will likely be difficult (Alexan- 
dei 1999 2005). The observed density profile in the inner 0.4 
pc falls as n(r) cx r~ 1A (Genzel et al. I 20031) . less steeply 
than predicted by the single mass BW76 solution. It is highly 
likely that additional dynamical processes beyond the simple 
picture considered here are at work, such as mass segregation 
(BW77). Furthermore, a large fraction of this observed dis- 
tribution is comprised of massive stars with lifespans shorter 
than either the NR or RR relaxation times (Fig. Q; these stars 
probably do not trace the relaxed population. 

We note that if the power-law profile were to continue to 
very small radii, the smallest radius r\ where a star is still sta- 
tistically expected to be found, r\ = N^ z ~ a \h ~6xl0~ 5 pc 
(for a= 1.4), would be well inside radius where RR depletes 
the cusp. If such a yet undetected population of tightly bound 
stars exists, it should exhibit an exponential central suppres- 
sion at ~ 10~ 3 pc. 

5.2.3. RR and dark mass in the Galactic Center 

It is likely that the Galactic MBH is surrounded by a dy- 
namically significant distribution of dark mass in the form of 
compact remnants, in particular SBHs that have sunk to the 
center by mass-segregation over the life time o f the Galaxy 
(Morris l 19931: M iralda-Escude & Gould l 20001 Ho pman & 
Alexander 2006; Freitag, Amaro-Seoane & Kalogera 2006). 
The existence of dark matter can be constrained by detecting 
orbital deviations from purely Keplerian motions (Mouawad 
et al. I 20041: Alexande r! 20051) . The co-existence of a dense 
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FIG. 8. — The density profile n(r) and the distribution of the semi-major 
axis n(a) in the GC for several values of \. 



cluster of SBHs with the dense stellar cusp may have in- 
teresting dynamical implications (e.g. Alexander & Livio 
2004). Efficient RR could deplete this dark mass component 
by rapidly draining the compact remnants into the MBH. For 
our assumed RR efficiency factor x = 1 me depletion will be 
significant only on scales < few x 10~ 3 pc (Fig. [3} with negli- 
gible dynamical implications for the observed stars. However, 
if \ > 10, the inner 0.01 pc will be depleted. Future high 
resolution observations of stellar orbits may therefore further 
constrain RR. 

RR may also affect the emission of gamma rays by annihi- 
lation of hypothetical dark matter (DM) particles near MBHs. 
It was shown that gravitational scattering between DM and 
stars le ads to a DM d ensity p rofile pdm oc r~ 3 / 2 (Gnedin & 
Primack 2004; Merritt 2004). This profile implies an annihi- 
lation rate which is proportional to log(r^/r; n ). In presence 
of RR, the inner cut-off of the cusp r m will be large as com- 
pared to the situation discussed in the literature where RR is 
absent, thus somewhat decreasing the annihilation rate. 

6. DISCUSSION AND SUMMARY 

RR is a coherent relaxation process that operates when sym- 
metries in the gravitational potential restrict the evolution of 
the orbits (R96; RI98; Tremaine 2005). As a result, the or- 
bits maintain their orientation over many orbital periods, and 
over that time the stars exert coherent mutual torques. These 
efficiently change the orbital angular momentum J (energy 
relaxation continues on the slow NR timescale). RR operates 
only under specific conditions. However, when these are met, 
it can be orders of magnitude more efficient than NR. Here 
we consider scalar RR, which affects both the direction and 
magnitude of J, in the regime around a MBH where the orbits 
are Keplerian (where both mass and GR precession are negli- 
gible), and vector RR, which affects the direction but not the 
magnitude of J, in the regime around a MBH where the po- 
tential is spherical. Scalar RR can deflect stars into the MBH, 
whereas vector RR can only change their orbital orientation. 

In this paper we studied the effect RR has on the energy 
DF of stars near a MBH. We explored the consequences for 
the disruption and capture of stars by the MBH, in particular 
through inspiral by GW emission, and used our results to in- 



terpret the properties of the observed dynamical components 
in the GC. The complex full 2+1 (£,J;t) problem was re- 
duced to an approximate 1+1 (£;t) Fokker-Planck equation, 
which we solved numerically. We also carried out Monte 
Carlo experiments to test some of our assumptions. 

Our results are as follows. (1) RR leads to the depletion 
of the high-energy end of the DF (stars tightly bound to the 
MBH), accompanied by an enhanced current of stars to high 
energies. The exact extent of the depletion and the effec- 
tive high-energy cutoff depend on the poorly determined ef- 
ficiency of RR. The currently available estimates of the RR 
efficiency factor indicate that x ~ 1. (2) We confirm the re- 
sult of RT96 that for % = 1, the direct tidal disruption rate 
is only modestly enhanced by RR. This is because scalar RR 
is most efficient close to the MBH where the orbits are Ke- 
plerian, whereas most of the tidally disrupted stars originate 
farther away from the MBH. A higher RR efficiency than as- 
sumed here could substantially increase the tidal disruption 
rate. (3) We show, in contrast, that the GW inspiral rate is 
dominated by RR dynamics, and is likely increased by almost 
an order of magnitude relative to the rate predicted assuming 
NR only. This is because stars undergoing GW inspiral orig- 
inate very close to the MBH. A higher RR efficiency than as- 
sumed here could substantially decrease the GW inspiral rate 
by completely depleting the tightly bound stars. (4) We apply 
our results to the GC and show that vector RR can naturally 
explain the inner cutoff of the ordered star disks at 0.04 pc and 
the transition to the randomized inner S-star cluster, as well as 
the randomized state of the old red giants in the inner ~ 1 pc. 
We also show that scalar RR is consistent with the presence 
of the disks, as it is too slow to disrupt them. If the S-stars 
were born in the disk on circular orbits, then scalar RR may 
explain may explain their present non-zero eccentricities. 

We note that there are additional processes that could be 
affected by RR, which we did not study here. Stars that un- 
dergo tidal capture and subsequent tidal heating ("squeezars", 
Alexander & Morris 2003) also originate close to the MBH, 
where RR dominates the dynamics. Squeezars could be di- 
rectly observed in the GC if RR substantially enhances their 
rates. A similar process is the tidal capture of a binary com- 
panion by an intermediate mass black hole in a young clus- 
ter and the subsequent Roche-lobe feeding that could power 
an ultra-luminous X-ray source (Hopman, Por tegies Z wart & 
Alexander 2004; Hopman & Portegies Zwart l"2005l Baum- 
gardt et al. I 20051) . The role of RR in tidal capture is still 
unclear, as there is no RR quenching by GR precession at the 
tidal capture radius that can stop the star from being rapidly 
destroyed. 

There are several limitations and uncertainties in our results 
that will have to be addressed by future studies. Our treatment 
of the problem in £-space incorporates the RR sink terms in 
approximate form only. The efficiency of RR is poorly de- 
termined as it has been calibrated based on a restricted set of 
iV-body simulations and its dependence on the parameters of 
the stellar system has not been investigated in full. This un- 
certainty is significant, as it could potentially reverse the sign 
of the effect of RR on the GW inspiral rate. While the rate is 
increased by an order of magnitude if % = 1, it is completely 
suppressed if % > 10 (Fig. 0. An important omission in 
our treatment is the assumption of a single mass stellar pop- 
ulation. A multi-mass population will induce mass segrega- 
tion. This will modify the DF directly (BW77), and likely 
also change the RR efficiency. The stellar DF also depends 
on processes that destroy one type of star while not affecting 



RESONANT RELAXATION 



11 



others , such as stellar collisions (e.g. Freitag & Benz l 200lt 
l"200l . 

Progress on these issues will likely require large scale nu- 
merical simulations. The effects of RR on the DF of stellar 
systems with MBHs have not yet been studied by A-body 
simulations. The recent full y self-consisten t A-body simu- 
lations of Baumgardt et al. f 2004a, 2004b|, Preto et al. ( 
2004) and Merritt & Szell f2005l) . indicate that such stud- 
ies are almost within reach of current hardware. Such simu- 
lations are particularly important for predicting GW inspiral 
rates, as these depend on a combination of complex mech- 
anisms, including mass-segregation and RR, which are hard 
to assess with (semi-) analytical methods (e.g. Baumgardt et 
al. I 20051) . Our analysis shows that it may be of consider- 
able importance to use a relat ivistic potential in such simu- 
lations (see e.g. Rauch 1999), since otherwise it is unclear 
whether RR is quenched, and a large number of stars may fall 
directly into the MBH without emitting an observable GW 
signal. Yet larger-scale modeling will require abandoning the 
exact A-body approach in favor of approximate methods such 



as Mo nte Carlo simulations with RR (Freitag & Benz l 200ll 
2002) or numerical 2+1 Fokker-Planck models. 

We note that throughout we assumed a spherically sym- 
metric DF. Deviations from spherical symmetry may affect 
both RR and the loss-cone structure in phase-space (see e.g. 
Magorrian & Tremaine [1999], who discuss the possibility 
of a loss 'wedge' in case of a triaxial DF). Typically this en- 
hances the event rates. 

Finally, we note that the observed dynamics of the GC pro- 
vide a promising empirical basis for calibrating and cross- 
checking theoretical and numerical studies of RR. 
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APPENDIX 

A. THE LOGARITHMIC FORM OF THE FOKKER PLANCK EQUATION 

Because of the large ratio of the tidal radius and the radius of influence, the natural way to integrate the Fokker-Planck equation 
is to divide the energy range into equal logarithmic intervals. For convenience we give here the equations in terms of the 
logarithmic distance variable z = ln(l + x). The Fokker-Planck equation Mil without sink terms is then written as 



where Q = T,f =1 Qi(z,r) with 



8t { > dz : 



g 1 (z,r)= 5 (z,r) 2 (e 2 -l)' 3 / 2 I 
Q 2 (z,r)=g(z) j*° dw{e w - l)~^^fl 



(Al) 



-e-'(e»-l)^ aflg( * ,r) 
dz 



dwe w g(w) , 



_ dg(z,r) I 1-1-3/2 w / \ 

— e — / aw(e — 1) ' e g[w) . 



The logarithmic expressions for the sink terms can be included directly from Eq. 1271 by replacing x 
the numerical integration is simple and can be easily done over many orders of magnitude in energy. 

B. THE NON-RESONANT LOSS-CONE SINK TERM 



(A2) 

1. Written this way, 



For convenience we reproduce here Eq. (14) of BW77, which interpolates between the empty and the full loss-cone regimes of 
the NR sink term, 



0.5g(x,r)x 5 / 2 



where 



aln(M./M*) [5.56 + ki(x D /4x)/q(x, r)] ' 
a = (M*/M.) 2 n h r\x D , q(x, r) = 1.6 In (6M./nM±) ax' 5/2 g(x, r) . 



(Bl) 



(B2) 
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